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We investigate the multi-wavelength emission of BzK selected star forming galaxies at z ~ 2 in the Great Ob- 
servatories Origins Deep Survey (GOODS) North region. Most (82%) of the sources are individually detected 
at 24/im in the Spitzer MIPS imaging, and one fourth (26%) in the VLA radio data. Significant detections 
of the individually undetected objects are obtained through stacking in the radio, sub mm and X-ray domains. 
The typical star forming galaxy with stellar mass ~ 10 ri M ( :>. at z = 2 is an Ultra-luminous Infrared Galaxy 
(ULLRG), with Lir ~ 1-2 x 10 1 L & and star formation rate SFR w 2OO-3OOM yr~ 1 , implying a comoving den- 
sity of ULIRGs at z = 2 at least 3 orders of magnitude above the local one. SFRs derived from the reddening 
corrected UV luminosities agree well, on average, with the longer wavelength estimates. The high 24/im de- 
tection rate suggests a relatively large duty cycle for the BzK star forming phase, consistently with the available 
independent measurements of the space density of passively evolving galaxies at z > 1 .4. If the IMF at z — 2 is 
similar to the local one, and in particular is not a top-heavy IMF, this suggests that a substantial fraction of the 
high mass tail ( >, 10 Mq) of the galaxy stellar mass function was completed by z ~ 1 .4. 

Subject headings: galaxies: evolution — galaxies: formation — cosmology: observations — galaxies: star- 
bursts — galaxies: high-redshift 



1. introduction 

The rate at which stars in massive galaxies were assembled 
is a crucial measurement for the characterization of galaxy 
formation. Hierarchical clustering models have traditionally 
favored galaxy formation at quiescent rates (e.g. Cole et al. 
2001). Formation of stars at high rates in massive galax- 
ies would qualitatively match monolithic formation scenarios 
(Eggen et al. 1962). Intense star formation at high redshift 
has been established for submm selected galaxies. These are 
relatively extreme objects, with space densities nearly two or- 
ders of magnitude smaller than local massive galaxies (e.g., 
Scott et al. 2002). Extending the census to less extreme star 
formation rate (SFR) levels is necessary, but obtaining reli- 
able SFR estimates for high-redshift galaxies is a challenging 
endeavor. The widespread presence of dust, absorbing the 
light of high mass stars and re-radiating it at longer wave- 
lengths, complicates the immediate use of the UV luminosity 
as a SFR indicator, as uncertain large corrections are required. 
A multi-wavelength approach, although observationally de- 
manding, can yield more robust measures of SFR. 

The major growth and assembly of galaxies' stellar mass is 
observed during the epoch between redshifts 1 < z < 3 (Dick- 
inson et al. 2003; Rudnick et al. 2003). Daddi et al. (2004b), 
on the basis of the highly complete spectroscopic database of 
the K20 survey (Mignoli et al. 2005), showed that a relatively 
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clean, efficient and complete (reddening independent) selec- 
tion of massive galaxies in the above redshift range can be ob- 
tained by selecting outliers in a (B — z) vs. (z — K) diagram. 
For z ~ 2 objects with K < 20, selected in this way, Daddi et 
al. (2004a;b; hereafter D04a and D04b) derived high masses 
and SFRs, suggesting they are the progenitors of local mas- 
sive spheroids caught during their phase of major assembling. 
The SFRs derived from the dust-corrected UV luminosities 
are however uncertain. 

In this letter, we have taken advantage of the BzK selec- 
tion to assemble a statistical sample of K < 20 (Vega) mas- 
sive star forming galaxy candidates at 1.4 < z < 2.5 in the 
Great Observatories Origins Deep Survey (GOODS) North 
field. Deep 24/im observations, obtained with MIPS on board 
of the Spitzer Space Telescope (SST), were used to study their 
rest-frame mid-IR emission properties. Deep X-ray, submm 
and radio data allow us to build a panchromatic view of their 
spectral energy distributions (SEDs), and to shed light on 
their nature and star formation activity. We discuss the im- 
plications for the assembly process of massive galaxies. We 
use a Salpeter initial mass function (IMF) from 0.1 to 100 
M , and a WMAP cosmology with S1 A , Q M = 0.73,0.27, and 
h = H [km s- 1 Mpc-']/100 = 0.71. 

2. BzK SELECTION OF z ~ 2 GALAXIES IN GOODS-NORTH 

The GOODS-North field has been observed in the /iT-band 
with the Flamingos camera at the Mayall 4-m NOAO tele- 
scope. About 4-7 hours integration were collected with aver- 
age 1 .2" seeing, over each of two contiguous pointings, reach- 
ing 5cr limits of K ~ 20.5 (Vega) for point sources. Calibra- 
tion was done using objects in common with 2MASS. Obser- 
vations in the B- and z- bands were obtained at the Subaru 
telescope with Suprime-Cam with < 1" seeing, and are de- 
scribed in Capak et al. (2004). Sources were selected over an 
area of 154 arcmin 2 to K < 20 (Vega), thus matching the limit 
where the BzK criterion is currently calibrated by the K20 sur- 
vey (D04b). The Capak et al. (2004) zeropoints were slightly 
adjusted, and a small color term was applied to the Z?-band 
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magnitudes because the Subaru B-band filter is redder than the 
one used in D04b. 169 z ~ 2 star forming galaxy candidates 
were selected having BzK = (z — K)ab — (B — z)ab > —0.2. 
We are explicitely excluding from the analysis the 13 objects 
with BzK < —0.2 and (z — K)ab > 2.5 that are candidate pas- 
sively evolving galaxies at z > 1 .4. However, the shallower 
depth of the B-band data (compared to those in D04b) pre- 
vents a clean separation of z > 1 .4 passive and star forming 
galaxies among sources with the reddest z — K colors. We 
might expect up to 10-15 genuine passively evolving z > 1.4 
galaxies among the BzK > —0.2 objects, counting the sources 
with no or low-significance B-band detection. 38/169 (22%) 
objects were discarded as likely AGN-dominated, because de- 
tected in the hard X-ray band (Alexander et al. 2003). The 
AGN contamination that we recover here is higher than in 
D04b, owing to the deeper X-ray data. The surface density of 
the 131 non X-ray bright BzK > -0.2 objects in GOODS-N 
is 0.85 ± 0.07 arcmin~ 2 (Poisson), consistent with that mea- 
sured in the K20 survey. On the basis of the K20 survey re- 
sults (D04b), we assume in the following that the 131 galax- 
ies with BzK > -0.2 are mostly at 1.4 < z < 2.5, with a 
relatively fiat redshift distribution in that interval and with 
a £ 10% contamination of objects at 1 < z < 1.4, and that 
our sample is complete for star forming galaxies with K < 20 
at 1.4 < z < 2.5. This is also supported by our photometric 
redshifts, and by the limited amount of spectroscopic redshift 
currently available for K < 20, BzK sources in GOODS-N. 

3. MULTI-A OBSERVATIONS OF BzK SELECTED GALAXIES 

Rest frame UV: The calibrations of D04b were used for 
estimating stellar masses, reddening and SFR of BzK se- 
lected galaxies from their observed properties in the opti- 
cal/IR. The average mass of our z — 2 object is 1 .0 x 10 1 l M & 
and the sample is complete above ~ 10 1 'M for 1 .4 < z < 2.5. 
The median B — z = 1.50 color translates to a reddening of 
E(B - V) = 0.40 for the Calzetti et al. (2000) law. This 
median reddening is similar to our estimate of the E(B — V) 
upper limit for the UV criteria of Steidel et al. (2004) for 
selecting z ~ 2 galaxies, and we estimate that roughly 50% 
of the z ~ 2 BzK star forming galaxies with K < 20 would 
be missed when selecting in the UV. We find an average 
SFR ~ 220 M Q yr _1 for our sample, based on the estimate 
of the reddening corrected 1500 A rest frame luminosity, with 
a median (average) correction factor of 46 (118). The ratio of 
the averages of the corrected and uncorrected SFRs is 42. 

MIPS 24\im: Deep 24fim observations of GOODS-North 
were obtained with SST MIPS, as a part of the GOODS 
Legacy Program (M. Dickinson et al., in preparation), for 
a total of 10.4 hours exposure time per sky pixel. Sources 
were detected in the 24^m data using SST IRAC data as prior 
positions, in order to improve source deblending (R. Chary 
et al., in preparation). We cross-correlated the 131 BzK se- 
lected galaxies to MIPS objects using a 2" search radius (ex- 
pected false matching rate of order of a few %). 107/131 
(82%) of the z ~ 2 galaxies are matched to a MIPS counter- 
part. The 24^m flux (/ 2 4) ranges from 300-500 [iJy for the 
brightest sources to the typical 3<r limits in the range of 15- 
30/Jy (depending on position) for undetected sources. The 
median /24 is 1 1 0/jJy, including non detected sources. The 
average is 127-124^Jy, depending on whether for the unde- 
tected sources we use the 3a upper limit or assign zero flux 
to them. A Kendall's r test indicates that /24 correlates with 
the Zf-band flux, with a 99% level of significance. Comparing 
to the full GOODS-N MIPS 24^m catalog, BzK star forming 



galaxies account for about 7-10% of 24/xm selected sources 
above 50 to 300 [iJy, with the fraction varying slightly with 
24/im limiting flux. These figures place lower limits to the 
fraction of z > 1.4 galaxies in 24/im selected galaxy sam- 
ples, which are consistent with the predictions of Chary & 
Elbaz (2001; hereafter CE01), i.e. 16%. At 1.4 < z < 2.5 
the 24/im observations probe rest-frame wavelengths of 7- 
10/im, the mid-IR region dominated by strong PAH emis- 
sion features and silicate absorption. In order to constrain 
the level of IR luminosity (L iR = L8-iooo>mX hence SFR, re- 
quired to reproduce the observed fan, levels, we used the mod- 
els from CE01, which provide LrR-dependent templates cal- 
ibrated from the local SEDs of IR luminous galaxies. For a 
given Lir, these models predict a factor of « 5 decrease of frA, 
from z = 1.4 to z = 2.5. Assuming a flat redshift distribution 
within 1.4 < z < 2.5 for the BzK > -0.2 selected galaxies, 
models with Ljr = 1 .7 x 10 12 L© are required to reproduce the 
125^Jy average flux level 1 . This corresponds to an average 
SFR~300M Q yr- 1 (Kennicutt et al. 1998). Using instead the 
empirical SED of Arp 220, calibrated by ISO observations of 
the PAH features region, one would expect a relative /24 peak 
around z ~ 2 and minima toward z = 1.4 and z = 2.5. This is 
due to strong 9.7/im silicate absorption, a feature which has 
been observed to be quite common at z <~ 2 in sources with 
Ljr ~ 1O 13 L (e.g., Yan et al. 2005). The / 24 expected from 
Arp 220 when averaged within 1.4 < z < 2.5 is about 30^Jy. 
For the Arp 220 IR luminosity of 1.5 x 1O 12 L0, this in turn 
implies Ljr ~6x 10 l2 L & for the typical BzK star forming 
galaxy at z = 2, or very large SFR ~ 1000 M Q yr _1 (this high 
value is disfavored by the other measurements, see below). 

Radio 20cm: 34/131 (26%) of the BzK > -0.2 galaxies are 
individually detected with S/N > 3 in deep VLA radio data 
(Richards 2000; reprocessed by G. Morrison et al., in prepa- 
ration), with fluxes at 1.4 GHz in the range of 20-160 /iJy 
(when excluding a radio galaxy with lmJy flux), and an aver- 
age of 37 [iJy. Stacking of the 94 undetected sources 2 , after 
discarding 3 objects close to unrelated radio sources, yields 
a S/N = 8 detection for an estimated average flux density of 
10 ±2 fiJy. When adding to this the individual detections 
we derive an average flux density of about 17yuJy for the full 
sample of 131 objects. Using a radio spectral index a = —0.8 
(/„ °< i/~ a ) we derive a luminosity of 3.6 x 10 23 W Hz -1 for 
an average < z >= 1 .9, which implies Lir ~ 1.2 x \Q 12 Lq and 
SFR ~ 210 Moyr" 1 (Yun et al. 2001; Kennicutt et al. 1998). 

SCUBA 850 [im: Only one of the 131 sources appears in 
the list of Pope et al. (2005) with a flux of 4mJy. In or- 
der to estimate the average 850/im flux of the undetected ob- 
jects we used the Borys et al. (2003) HDF-N SCUBA su- 
permap. Some 97 MlPS-detected objects were stacked 3 , af- 
ter excluding those separated by less than a SCUBA beam 
halfwidth (7") from any known SCUBA detection. As the 

1 Note that we are studying the 24^tm properties of a near-IR selected 
sample of galaxies, not a 24/im selected sample. Its redshift distribution is 
expected to be flat within 1.4 < z < 2.5, independent of the behavior of the 
24/im flux with redshift. As we are detecting 82% of the objects, the typical 
flux that we detect corresponds to the source at typical redshift. 

2 To stack the radio data, sub images were extracted at the locations of the 
non-detected sources. These were corrected for VLA primary beam attenua- 
tion and combined using a weighted average. The integrated flux density and 
error were computed with the AIPS task JMFIT, which modeled the stacked 
emission using an elliptical Gaussian. 

3 We are using only MlPS-detected objects to maximize the signal to noise 
by avoiding passive galaxy contaminants, while the possible star forming 
galaxy contaminants at 1 < z < 1.4 would not alter significantly the stack- 
ing because of the mild dependence between 850^tm flux with redshift. 



Multi-A observations of z — 2 galaxies in GOODS-N 
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FIG. 1. — The average multi-wavelength emission of 131 BzK selected 
star forming galaxies at 1.4 < z < 2.5. The model shown here (from CE01) 
has Lir = 1.7 x 1O 12 Z,0 and is redshifted to z = 1.9. Measurements in SST 
IRAC ch2 (4.5/tm, AB = 20.78 mag) and ch4 (8.0/im, AB = 20.99 mag) were 
obtained from the data of M. Dickinson et al., in preparation. The horizontal 
bar on the right shows the 2-8 keV emission that would be expected from a 
typical AGN having the mid-IR flux observed for our objects. 



depth of the SCUBA data vary considerably over GOODS- 
N, we used weights in the stacking based on the local noise. 
The average 850/im signal is 1.0±0.2mJy. We checked that 
this 5-sigma detection is robust using Monte-Carlo simula- 
tions. Scaling down the result by 18% to account for the MIPS 
undetected objects, the average 850/im flux corresponds to 
Lir ~ 1.0 x 1O 12 L , both for the CE01 templates and when 
using Arp 220, or to SFR ~ 170 M yr _1 . At 450/mi we ob- 
tain a 3cr upper limit of 6mJy, consistent with the models. 

Chandra soft and hard X-rays: Is the detected IR emission 
AGN or starburst powered ? At z ~ 2 the Chandra hard 2- 
8 keV band is sensitive to highly penetrating rest-frame X- 
rays up to about 20 keV, which can escape large column den- 
sities up to about 10 24 cm~ 2 . The lack of hard X-ray detec- 
tion in the Chandra 2 Ms exposure should ensure that bright 
AGNs are generally not included in our sample. The typi- 
cal upper limits to the 2-8 keV flux of order < 3 x 10~ 16 
erg s cm~ 2 coupled to the typical 24/im fluxes measured, 
imply lower limits to the 24/im to 2-8 keV flux ratios, for 
single objects, about one full order of magnitude larger than 
expected from AGNs (see e.g. Rigby et al 2004). This sug- 
gests that the 24/im emission in our sample is powered by 
star formation in most cases. While we have discarded all 
direct hard X-ray detections from our sample, 6/131 (4.5%) 
of the objects are detected in the soft-band only, with fluxes 



?i 10 erg s cm (Alexander et al. 2003). Four of these 
are detected also in the radio, and one is the only SCUBA 
detection. These soft sources are most likely the most ex- 
treme starbursts in the BzK > —0.2 sample, similar to object 
ID#5 described in D04a. We stacked the X-ray undetected 
sources, considering only the GOODS-N region within 6 ar- 
cmin of the 2 Ms Chandra data aim point where the sensitivity 
is the highest, and avoiding all objects closer than 3 Chandra 
PSF from known X-ray sources. The stacked images result 



in a 8.5(7 detection in the soft 0.5-2 keV band, for a flux of 
1.0 x 10~ 17 erg s cm -2 . The non detection in the hard 2- 
8 keV band constrains the photon spectral index T > 1 .0 at 
the 3-sigma level, and implies an average 24/im to 2-8 keV 
flux ratio over two orders of magnitude above what is typical 
for AGN at z ~ 2 (Rigby et al. 2004). If AGN are present they 
would have to be heavily Compton thick. Using T = 2.0 as ap- 
propriate for starbursts implies a rest-frame luminosity in the 
2-10 keV range of 3.4 x 10 41 ergs" 1 for <z>= 1.9. Adding 
back the individual soft X-ray detections would increase the 
2-10 keV luminosity to about 5 x 10 41 erg s _1 . This is a fac- 
tor of two lower than estimated by D04b for K20 BzK sources. 
Using the Ranalli et al. (2003) calibration this translates into 
an average SFR ~ 100 .l/.yr '. or SFR ~ 500 Msyr" 1 if 
using the Persic et al. (2004) calibration. The large differ- 
ence is due to the relative expected importance of high ver- 
sus low mass X-ray binaries. The latter calibration appears 
more appropriate when dealing with IR luminous sources with 
Lir > 1O 12 L (Persic et al. 2004). 

4. DISCUSSION 

The multi-wavelength SED (Fig.Q of BzK selected z — 2 
star forming galaxies consistently indicates an average Lir ~ 
l-2xlO 12 L and SFR ~ 200-300 M yr-\ supporting the 
earlier claims of D04a;b and implying that the local correla- 
tions (e.g., the mid and far-IR to radio correlations) hold, for 
the average BzK galaxy, also at z ~ 2. While the X-ray based 
estimate appears the least accurate, due mainly to the large 
uncertainties in its calibration, the X-ray properties clearly 
support that the mid to far-IR emission of these sources is in- 
deed dominated by vigorous star formation and not by nuclear 
activity. The inferred average Lir £ 1O 12 L implies that the 
typical BzK selected star forming galaxy is an Ultra Luminous 
IR Galaxy (ULIRG). Morphology from HST ACS indeed sug- 
gests that, in many cases, these z = 2 galaxies are assembling 
through merging (D04a), similarly to local ULIRGs. Using 
the volume at 1.4 < z < 2.5 (5.7 x 10 5 Mpc 3 ), we can put a 
lower limit to the spatial density of ULIRGs at z ~ 2 of about 
1-2 x 10~ 4 Mpc -3 . This is about 3 orders of magnitude higher 
than at z ~ 0. 1 (Sanders et al. 2003), and a factor of 2-3 higher 
than at z = 1 (Le Floc'h et al. 2005). While ULIRGs are ex- 
ceptional objects for the local universe, they appear to be the 
norm among massive z ~ 2 star forming galaxies. 

From the X-ray we also derive a fairly low X-ray to op- 
tical flux ratio of log (fo. 5-2 kev/fR) ~ —1-7, which suggests 
that also the rest frame UV emission of these z = 2 galax- 
ies is dominated by the emission of stars. The average SFR 
inferred from the dust-corrected UV luminosity agrees well 
with the longer wavelength estimates. Fig.|2](top) shows the 
ratio of 24/im to B-band flux densities as a function of the 
rest-frame UV galaxy colors. Apart from K-correction ef- 
fects, this observed mid-IR to UV luminosity ratio measures 
the ratio of dust-extinguished to relatively unobscured SFR. 
The Kendall's r test detects a correlation at > 99.9 confi- 
dence level. A similar trend is found when using the radio 
flux instead of f^. This confirms that, in general, the red UV 
continua of BzK selected z = 2 star forming galaxies are in- 
deed due to dust reddening. The diagonal line plotted in Fig. 13 
(top) shows the expected scaling in the case that the UV and 
mid-IR trace the same amount of SFR. The average ratio and 
its dependence with the UV color are in line with what is ex- 
pected. Sources with higher 24/im flux, however, tend to have 
systematically larger ratios, and vice-versa. This can be more 
clearly seen in Fig. [2] (bottom), where we plot the 24/xm to 
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FIG. 2.— (Top pane/). The ratio of 24pm to B-band flux for K < 20 Bz£ 
selected z = 2 star forming galaxies in GOODS-North as a function of B — z 
color (i.e. dust extinction at 1500A, as inferred using a Calzetti et al. 2000 
law). The diagonal line shows the expected ratio for the average z = 1.9 
in the case that the UV and mid-IR trace the same amount of star formation. 
Solid and empty symbols are for sources with /24 above or below the median, 
respectively. (Bottom panel). The dust extinction corrected 24^m to B-band 
flux ratio, as a function of the 24pm flux. Solid and empty symbols are for 
sources with UV reddening below or above the median, respectively. In both 
panels, values on axis labeled with Lir or SFR are for an average z = 1 .9 and 
were computed using the CE01 model reproducing the average SED (FigEJ. 



B-band flux ratio (with the B-band flux here being corrected 
for dust reddening) as a function of the 24/im flux: a clear 
trend with 24/im flux is present. The expected K-correction 
term is small, and this plot suggests that sources with higher 
Ltr have progressively larger fractions of the SFR that cannot 
be recovered from the UV luminosity even after correcting 
for dust reddening. However, such a strong trend is not re- 
produced using the radio data for the radio detected objects. 
Fig. |3 (bottom) also shows that, at fixed /24, the sources with 
bluer UV continua tend to have higher corrected 24/xm to B- 
band flux ratio by a factor of 2 on average. 

Vigorous starbursts are present within the BzK selected star 
forming galaxies at z = 2, which have also fairly large typical 
stellar masses of - lO n M for # < 20 (Vega). These masses 
would grow even larger as a result of star formation, depend- 
ing on the duty cycle of the star formation event. The high de- 
tection rate (82%) at 24/im supports the possibility that high 



SFRs among these K < 20 sources are sustained for a sub- 
stantially long amount of time during 1.4 < z < 2.5. When 
limiting to a common stellar mass threshold, e.g. > 10 n M Q , 
the space densities of passively evolving galaxies is compa- 
rable to that of vigorous starbursts within 1 .4 < z < 2, and 
perhaps much smaller at 2 < z < 2.5 (Daddi et al. 2005; Kong 
et al. 2005; McCarthy et al. 2004). This would suggest that 
the average duty cycle is at least 50%, and likely more. The 
Universe ages from 2.6 Gyr to 4.6 Gyr between z — 2.5 and 
1 .4, and about 1 Gyr is still available on average per galaxy 
before z = 1 .4, implying a continuation of the star formation 
event for order 0.5 Gyr or more, on average, and that typically 
these objects have been active for a similar or larger amount 
of time before observations. This is also consistent with the 
typical age of the present SF event inferred from the opti- 
cal/IR SED fitting with constant star formation rate models, 
which is mainly based on the strength of the observed Balmer 
break (about 0.7 Gyr; D04a; see also Shapley et al. 2005). 
Therefore, by z — 1.4 the typical mass of these galaxies will 
have roughly doubled, reaching — 2 x 10 M©, on average. 
From their observed space densities (~ 2.3 x 10~ 4 Mpc~ 3 ) 
and SFRs, we infer that the integrated stellar mass density 
formed within BzK > —0.2 galaxies with K < 20 in the 2 Gyr 
time within 1.4 < z < 2.5 is ~ 10 8 M Q Mpc~ 3 , independent 
of the duty cycle. This is only — 20% of the local total stel- 
lar mass density, but is comparable to the local stellar mass 
density for objects with stellar mass > 2 x 10 M© (Cole et 
al. 2001). If the IMF at z = 2 is similar to the local one, 
and in particular is not a top-heavy IMF, this suggests that 
by z ~ 1.4 the assembly of the high-mass tail (> lO n M ) 
of present day's galaxies mass function was, in a significant 
part, completed. This is also supported by the measurements, 
quoted above, of a comparable space densities of old and pas- 
sive galaxies with similar masses to the star-forming BzKs 
already existing at 1.4 < z < 2, and with the evidences (Pa- 
povich et al. 2005) of strongly declining specific SFRs for 
z < 1.4 massive galaxies. 
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